Parity in the CMB: Space Oddity 
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We search for a direction in the sky that exhibits parity symmetry under reflections through a 
plane. We use the natural estimator, which compares the power in even and odd £ + m multipoles, 
and apply minimal blind masking of outliers to the ILC map in order to avoid large errors in the 
reconstruction of multipoles. The multipoles of the cut sky are calculated both naively and by 
using the covariance inversion method and we estimate the significance of our results using ACDM 



simulations. Focusing on low multipoles, 2 < £ < 



with 



5, 6 or even 7, we find two 



perpendicular directions of even and odd parity in the map. While the even parity direction does 
not appear significant, the odd direction is quite significant - at least a 3.6a effect. 



I. INTRODUCTION 

Since the release of the COBE data and especially af- 
ter the release of the WMAP data many groups argued 
that there are unexpected features at the largest possi- 
ble scales (see p], and references within), typically at 
2 < £ < 5, (for a partial list, see |3|— ISI] ) - A real anomaly 
at the largest scale, if exists, could be explained by pre- 
inflationary physics. Thus having a better understand- 
ing of these potential anomalies might tell us something 
about the universe before inflation. 

Cosmic variance and foregrounds make it harder to 
estimate the significance of the various claims in the lit- 
erature. Moreover, a standard argument against most 
of these potential anomalies is the lack of motivation to 
look for them to begin with. Or as has been often re- 
marked (see e.g. [1]), when applying enough statistical 
estimators to a given data set with random behavior, 
"anomalies" are bound to appear. Estimating their sig- 
nificance, one tends to ignore this fact and focus on the 
anomalous finding alone. 

With this concern in mind we focus in this paper on a 
simple question that needs no justification to be asked: 
does the CMB at the largest scales behave as expected 
with respect to a mirror parity transformation? This is 
not the first time this question was asked. The authors 
of [E| studied the parity transformation of the CMB with 
respect to the direction of the alignment of the £ = 2,3 
multipoles Recently, [?J studied the parity of the 
CMB under reflections through the galactic plane. Here 
we study the parity symmetry of the CMB data with re- 
spect to all directions in the sky and check for anomalous 
parity directions. 

The paper is organized as follows. In section [III we de- 
fine a natural estimator for parity symmetry with respect 
to reflections through a plane normal to some direction n 
and present the masking scheme we use. In section HIJ we 
generate parity maps for the full and cut WMAP seven- 
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year Internal Linear Combination (ILC) map that show 
the parity score as a function of n. We find a parity-even 
direction that is related to the alignment of the £ = 2,3 
multipoles and is not significant. More interestingly there 
is a significant (about 3.6a) parity-odd direction. In sec- 
tion |IV] we conclude. 



II. METHOD 

A. Full Sky Parity Estimator 

We would like to check the CMB temperature fluctua- 
tions map for parity symmetry with respect to reflections 
through a plane, i.e. r r — 2 (r • n) n, where n is the 
normal to the plane. Using the full sky it is most conve- 
nient to work with the spherical harmonics Yg m (n) which 
transform as (— l)^+ m under this operator. Hence a nat- 
ural parity estimator is 



= E E {-if +m ^ 



>(n)r 



£=2 m= 



Ct 



(i) 



where n is the z-axis used in the harmonic expansion and 

The more S(n) is positive (negative) the more the 
temperature map is parity-even (odd) in the n direc- 
tion. Since the ensemble average does not vanish, (S) = 
tax — 1, it is convenient to redefine the parity estimator 
as 



S(n) = S(h) - 



!)■ 



(2) 



We shall use this estimator in the next section. 

This simple parity estimator is useful when dealing 
with full sky data. When working with actual data which 
suffers from foreground noise, one needs to mask some 
parts of the sky. The generalization of this parity esti- 
mator to this case is the subject of the rest of this section. 
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FIG. 1: The squared smoothed ILC map. The prominent 
regions lie close to the galactic plane, or inside it. 



FIG. 2: Masks that cover the most powerful pixels of the 
squared, smoothed ILC map. The masks cover 2.5% (black), 
5% (red), 7.5% (orange) and 10% (yellow) of the sky. The 
KQ75 galactic mask is also shown (in light gray). 



B. Masking Scheme 

Naturally, we would like the results to be insensitive 
to foregrounds. The cleanest approach would be to use 
the full WMAP KQ75 or KQ85 masks. Unfortunately, 
these masks cover roughly 30% and 20% of the sky, re- 
spectively. Such extremely large masks induce significant 
errors in the calculations of the multipole coefficients re- 
gardless of the reconstruction method used to calculate 
them and thus render the parity estimator unreliable. 

In an attempt to overcome this issue, we would like 
to define a mask that is as small as possible (and so 
it can be used to calculate multipole coefficients) but is 
still efficient in removing local outlying regions and so 
enables to test large scale physics. We emphasize that 
the mask defined below is by no means an alternative to 
the standard masks used to extract small scale properties 
of the CMB. 

To pinpoint the regions that induce noise at large scales 
we examine the squared temperature map after smooth- 
ing it by 20°, as shown for the ILC in Fig. [1] This is mo- 
tivated by [1[, where this map of anisotropy power was 
used to manually define localized regions and check the 
sensitivity of the alignment of £ = 2,3 to their removal. 
Here we simply automate this procedure to choose the 
most outlying local regions in the map. Not surprisingly, 
it is readily seen that the most intense areas in the map 
lie in small patches around the galactic plane. One op- 
tion would be to mask out only the intense regions that 
are inside the KQ75 mask. This option, however, is prob- 
lematic when estimating the significance of the result vs. 
random simulations of ACDM. The reason is that such an 
estimation requires an automatic procedure for defining 
the masked regions. But the overlap between the most 
intense regions of a random map and the KQ75 mask is 
very different than the overlap between the most intense 
regions of the ILC and the KQ75 mask. 

An alternative option, which we do follow, is to mask 
out a fixed total area of deviating pixels. For our anal- 
ysis we examine the masking of the most intense areas 
from 1% to 20% of the sky. For the ILC map, masks 



covering 2.5%, 5%, 7.5% and 10% of the sky are shown in 
Fig. [2j together with the galactic KQ75 mask, for refer- 
ence. Most of these masked regions lie inside or near the 
KQ75 galactic mask. 1 To quantify this, we calculate for a 
given masking area, A, the area that lies within the KQ75 
mask, A*, and compare this to random ACDM simula- 
tions. We find that for 5% < A < 15%, only -0.6% of 
the simulations have a greater A* than that of the ILC 
map. This provides some justification for our choice of 
masking these local regions around the galactic plane in 
the case of current WMAP data. 

Generally speaking there are two ways to calculate the 
multipole coefficients a£ m on a masked sky (both break 
down if the mask is too large). The naive method is 
to zero the region inside the mask and take a spherical 
harmonic transform. A more precise method is the 'co- 
variance inversion method' used in the analyses of [9l4ll) 
which we discuss next. 

The relation between the temperature map and the 
spherical harmonic coefficients is 

x = Ya + n , (3) 

where x is a vector of the temperature in the masked 
sky, Y is a matrix of the spherical harmonics evaluated 
on each direction (Y^- = 7^. mj (r^)), a is the coefficients 
vector and n represents the noise in the data. Assuming 
Gaussianity and statistical isotropy and ignoring the de- 
tector noise, the familiar solution for the reconstructed 
coefficients a is 

a= (Y+C^Y)" 1 Y+C^x. (4) 

The square matrix C is constructed from a given a power 



These regions overlap with the regions identified in Q as the 
ones containing most of the large scale angular power in the ILC 
behind the KQ75 galactic mask. 
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spectrum Ct as 

c ^= E 2 -^ p dii-ii)Ci, (5) 

^=^max + l 

where Pi is the Legendre polynomial. The sum in (j5j) 
runs over those multipoles that are not included in the 
coefficients vector. In this manner, this method compen- 
sates for the contribution of higher-^ multipoles that leak 
into the reconstructed coefficients since the spherical har- 
monics are no longer orthogonal on the cut sky. For the 
analysis in this work we use the seven year WMAP power 
spectrum and take L = 196, roughly corresponding to a 
scale of 1°. 



C. Cut Sky Parity Estimator 

Using the masking scheme discussed above, we replace 
our full sky estimator S(h) with a corresponding cut sky 
version 

S(n,A) = 

£ E (-l)^ '^y - (4>ax " 1) , (6) 
£—2 m——£ ^ ' 

where a^ m (n, A) and Ci(A) are reconstructed after mask- 
ing the area A of the most intense pixels in the map. 

If there is a direction in the sky, no, which is parity- 
even (odd) in an anomalous way then we expect S^rio, A) 
to increase (decrease) significantly as we increase A, since 
more and more galactic noise is removed. Eventually, as 
we keep on increasing A, the reconstruction noise grows 
and at some point it is bound to dominate and lower 
(raise) S(ho^A). In the next section we check if there 
are such directions in the WMAP data. We will mainly 
use the covariance inversion method. However, since the 
naive method is much faster to compute and does not 
rely on a particular power spectrum, we shall use it as 
well, for comparison. 



III. RESULTS 

We begin by calculating the parity score on the full sky 
WMAP seven-year ILC map degraded to HEALPix reso- 
lution A^ide = 8 after removing the Doppler-Quadrupole 
contribution [12|. The result for £ max = 5 (motivated by 
SB) is shown in Fig. [3] As can be seen, the maximum 
lies at (Z,6) = (260°, 48°) which is exactly the direction 
found in [4[ of the renowned alignment of the quadrupole 
and octupole. Since the alignment of the multipoles is di- 
rectly related to their planarity (corresponding to large 
m = ±£ modes), this can be seen as a special case of an 
even parity symmetry (large even £ + m modes) Q . On 
the full sky, this direction does not seem significant. 




FIG. 3: The parity estimator calculated on the full sky 7-year 
WMAP ILC map, with ^ ma x = 5. The maximum lies in the 
direction of the "axis-of-evil" . The score map is at HEALPix 
A S ide = 64 resolution, corresponding to 1°. 



We also notice a minimum in the map at (Z,6) = 
(266°,— 19°), signifying odd parity. However, this mini- 
mum as well does not appear to be significant on the full 
sky map. 

As discussed in [1], the alignment of the quadrupole 
and octupole is extremely sensitive to the masking of 
small patches near the galactic plane and even with small 
percentages of the sky removed, the alignment is de- 
stroyed. To check the stability of the parity results, we 
use the masking scheme described above and check the 
behavior of our estimator on the cut sky. When degrad- 
ing the cut sky map to = 8, we follow the prescrip- 
tion of whereby each downgraded pixel is calculated 
by averaging only the corresponding unmasked original 
pixels in the higher resolution, unless more than half of 
them are inside the higher resolution mask, in which case 
the downgraded pixel is masked as well. 

The results using two representative masks covering 
the most powerful 10% and 12% of the sky (both of 
which, as any masking above 2.5%, eliminate the align- 
ment) are shown in Fig. [H Comparing the results to 
the full sky map (Fig. [3|), we can see that the loca- 
tion of the maximum has shifted by almost 40°, to 
(Z,6) = (198°, 55°). It now stands out relative to the 
other positive peaks, but it does not appear significant. 

The result for the minimum of the map (which signifies 
maximal odd parity) is more compelling. Its location has 
shifted by no more than a few degrees when applying the 
masks and remains stable under further cuts (including 
the full KQ85 or KQ75 masks). Furthermore, it now 
appears much more significant, as we discuss below. 

It is worthwhile to mention that using the cut sky, 
the resulting even and odd extrema lie at a distance of 
90° ± 1° from each other. 

The final step in our analysis is the estimation of the 
significance of our findings. We have seen that on the 
masked ILC map, we get distinct even and odd parity 
directions. Noticeably, these directions in the score map, 
apart from having a high absolute value, also stand out 
relative to the rest of the map. To quantify this, we take 
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FIG. 4: The parity estimator calculated on the ILC map with part of the sky masked as described in the text. The masks 
used in panels (a) and (b) cover 10% and 12% of the sky, respectively. The multipole coefficients are reconstructed using the 
covariance inversion method. 



the maximum and minimum of the score map, 

(A) = max S (n, A) , S- (A) = min S (n, A) (7) 



and define 



S±(A) = 



S ± (A)-n(A) 



a(A) 



(8) 



where ji{A) and cr(A) are the mean and standard devia- 
tion of the score map, respectively. We can now compare 
these standardized scores to the values of S±(A) calcu- 
lated for random simulations. 

We generate random maps using the seven year 
WMAP power spectrum (at a resolution corresponding 
to the 1° resolution of the ILC) and degrade them to 
Aside = 8 before calculating the score. In order to com- 
pare with the results on the masked ILC, we apply the 
same masking procedure described above to the random 
maps. 

We find that the even parity direction (7,6) = 
(198°, 55°) is the maximum of the score map for mask- 
ing areas larger than 7%. However, it doe s not appear 
significant. This can be seen from Fig |5(a)| As discussed 
above if there is a direction in the sky which is parity-even 
in a significant way then S+(A) is expected to increase 
considerably and consistently as we increase A and even- 
tually decrease as we keep on increasing A. This is not 
the case in Fig 5(a) Indeed, if we compare the typical 
value of S+(A) of ~ 3 to random simulations, 12% get a 
higher score. 

The odd parity direction (l,b) = (26 4°, -1 8°), how- 



ever, is significant. As can be seen in Fig. 5(b) , when us- 
ing the covariance inversion method we see that S-(A) 
behaves exactly as expected from an anomalous direc- 
tion: it increases as we increase A until A is too large 
and the reconstruction noise takes over. Compared to 
random simulations, the significance of the peak score 
at A = 12% is 3.6<r, as 0.03% of 500,000 randoms get a 
higher score. 

When using the naive reconstruction method, we see 
that the score S-(A) follows the same behavior and 



reaches roughly the same maximum value as the covari- 
ance inversion method score. As expected, it peaks at 
lower masking areas. The reason for this is that the more 
robust covariance reconstruction compensates for the re- 
moved pixels using the given ACDM correlations, while 
the naive method completely ignores them. Therefore, in 
order to fully remove the influence of noisy regions, more 
masking is required when using the covariance inversion 
method. The discrepancy between the two methods sup- 
ports the choice of the robust method for estimating the 
significance of the results, which we have done using the 
small masks described above. 

It is important to note that the odd parity direction 
remains significant not only for the choice of ^ max = 5, 
but for 6 and 7 as well. In fact, for ^ max = 6 it reaches 
as high as 4.3<r. For larger values of ^ max the significance 
starts to drop. However, this is not surprising since the 
errors in the reconstruction method grow quickly with 
i [9l-[ll| . Moreover, as discussed below, if the origin of 
the anomalous parity behavior is pre-inflationary physics, 
then it is expected to drop as we increase t. 

One might worry that our masking scheme, the recon- 
struction methods or even our choice of degrading the 
maps to HEALPix resolution Ag^e = 8 before calculat- 
ing the score induce bias on the parity score that renders 
the parity-odd direction so significant. To check the ef- 
fect on the score and verify that it should not change our 
conclusions, we plot in Fig. [6] histograms for the score 
S- based on 10 4 random simulations. On each of these 
random maps three versions of the score were calculated: 
without masking and with Ag^e = 8, with masking of 
12% and N^&q = 8 and finally without masking but with 
A^side = 16. As can be seen, all three histograms are vir- 
tually identical. We therefore conclude that the masking, 
as well as the choice of HEALPix resolution, do not in- 
troduce meaningful bias to the score, and do not alter 
our result. In addition, very recently [14] discussed the 
covariance inversion method for reconstructing the spher- 
ical multipoles and made different choices in smoothing 
the initial sky map and in the parameters for degrading 
the map and cutting off the higher multipoles. We veri- 
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FIG. 5: The scores §+ (a) and S- (b) as a function of the masking area, A, using the covariance inversion method. The naive 
result is also plotted, for comparison. 
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FIG. 6: Histograms for the score S- calculated on 10 4 random 
realizations, without masking and with iV s ide = 8 (blue), with 
masking of 12% and iV s id e = 8 (green) and without masking 
and with iV s ide = 16 (red). The dashed vertical line is at the 
score value of 4.3, attained by ILC map. 



fled that the effect of such different choices on our results 
presented here is negligible as well. 

Finally, to support the possibility of a cosmological 
origin for our findings, we check our cut sky estimator 
directly on the WMAP V and W frequency band maps 
as well. Since in these maps the galactic foregrounds 
are not removed as in the ILC map, we cannot use the 
masking procedure described above. Instead, we use the 
KQ85 galactic mask (covering ~ 20% of the sky) for the 
reconstruction of harmonic coefficients, potentially intro- 
ducing significant reconstruction errors. Nevertheless, as 



is evident in Fig. [7J both bands show very similar odd 
and even parity signals, at the same locations. 



IV. CONCLUSIONS AND DISCUSSION 

We have found strong evidence for odd parity in the 
WMAP data. In addition, there appears to be a parity- 
even direction at 90° from it that does not appear signif- 
icant. These results were achieved using a novel masking 
scheme, that, as far as we can see, is the cleanest way to 
calculate large scale multipole coefficients from WMAP 
data. This masking scheme is a compromise between 
reducing remaining outliers in the low foreground ILC 
full sky map and avoiding reconstruction noise and is far 
from being ideal. It keeps most of the regions masked by 
the KQ75 mask and it ignores some regions outside the 
KQ75 mask. 

Planck is expected to see better than WMAP though 
the Galactic noise. Thus the significance of the parity- 
odd and even directions we found is likely to change. 
Since we are already reporting at least a 3.6a effect, if 
the significance increases further, we believe that this will 
be a real challenge to ACDM, one that calls for either a 
systematic or cosmological explanation. 

One possible cosmological explanation, that in fact mo- 
tivated us to look for non-trivial parity effects in the 
CMB, is a moving pre-inflationary particle (PIP). Mo- 
tivated by [HI the cosmological imprints of a static PIP 
were studied in [l6|, [l^ • The most distinct CMB imprint 
of a stationary particle was found to be giant concentric 
rings. Such rings were found in [l8| with fairly high sig- 
nificance (about 3.1<r). There is, however, no reason why 
a PIP should be static in the CMB frame. This moti- 
vated us to study the cosmological imprints induce d by 
the velocity of the PIP. As will be reported elsewhere [19j , 
the most general effect induced by the PIP's velocity is 
a non-trivial parity structure in the CMB. 
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FIG. 7: The parity estimator calculated with the covariance inversion method on the V and W frequency band maps with the 
KQ85 galactic mask applied is shown in panels (a) and (b), respectively. We also mark with a black dot the location of the 
ecliptic pole, at (/, b) = (276°, —30°), where the parity score is neutral. 



Another possibility is a pre-inflationary string. As far 
as we can see at the moment such an object explains 
naturally the parity-even direction. However, it is hard to 
explain the more significant parity-odd direction (and the 
giant rings of [18]) via a pre-inflationary string. It should 
be interesting to see if other pre-inflationary scenarios, 
such as bubble collisions (see [20| and references within), 
could induce the parity effects reported here. 



crowave Background Data Analysis (LAMBDA) (2lj and 
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